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Using Monte Carlo techniques, we show that the two-dimensional anisotropic Heisenberg model 
reproducing nicely inelastic neutron scattering measurements on La5CagCu24 04i (Matsuda et al. 
[Phys. Rev. B 68, 060406(R) (2003)]) seems to be insufficient to describe correctly measurements 
on thermodynamic quantities like the magnetization or the susceptibility. Possible reasons for the 
discrepancy are suggested. 
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Recently, Matsuda et alX reported results of ca- 
reful inelastic neutron scattering experiments on Las 
CagCu2404i. From standard spinwave analysis, it is con- 
cluded that the experimental data can be well reproduced 
by a two-dimensional anisotropic spin-1/2 Heisenberg 
model with antiferromagnetic nearest-neighbor interac- 
tions, J c i (= 0.2meV), and ferromagnetic next-nearest- 
neighbor couplings, J C 2 (= — 0.18meV), between the Cu 
ions in the CuC>2 chains, as well as antiferromagnetic 
couplings Jad (= 0.681meV) and J ac i {— 0.5J ac i = 
0.3405meV) between nearest and next-nearest neighbor- 
ing Cu ions in adjacent chains, respectively (see Fig. 1 
in Ref. Q). The uniaxial anisotropy of the spins along 
the b axis is written in the form of a single-ion inter- 
action D (= — 0.211meV), summing over the anisotropy 
contributions of the different couplings. 

This model reproduces the measured dispersion rela- 
tions along the Cu02 chains, i.e. along the c axis, and 
perpendicular to these chains, along the a axis. It is sup- 
posed to describe, at low temperatures, the long-range 
ferromagnetic order in the chains and an antiferromag- 
netic order perpendicular to the chains. 

We performed Monte Carlo simulations on the model 
in its classical limit with unit vectors as the spins, 
where the uniaxial anisotropy along the b axis is de- 
scribed either by a single-ion anisotropy)- 1 DJ2(Sj) 2 , or 
by anisotropic exchange interactions distributing the to- 
tal anisotropy in various ways among the different cou- 
plings. We used exactly the interaction strengths deter- 
mined in Ref. QJ- To compare the simulational data to 
previous measurements^ on the magnetization and spe- 
cific heat of La5CagCu24C>4i, we included external fields 
H oc along the b and c axes, and we recorded especially 
sublattice magnetizations, the total magnetization, the 
susceptibility, and the specific heat. We took care that 
reliable equilibrium data were obtained, using at least 
5 x 10 5 Monte Carlo steps per spin in each run. We stud- 
ied finite-size effects, simulating square lattices with the 
linear dimension ranging from 10 to 100. The square lat- 
tice is constructed in such a way that its principal axes 
are given by the two directions along which J ac \ acts 
between a reference site and the two sites in one of the 
adjacent chains. Thence the relatively strong coupling 
J QC i connects nearest-neighbor spins on that lattice, and 



the Cu02 chains run diagonally through the lattice. This 
geometry shows most directly the relation of the model 
by Matsuda et al. to much studied anisotropic Heisen- 
berg models with nearest-neighbor interactions on square 
lattices. 

Our main findings may be summarized as follows: 

(i) For vanishing field and at low temperatures, one 
obtains, indeed, antiferromagnetic order, i.e. the model 
describes a ferromagnetic ordering of the spins along the 
Cu02 chains and an antiferromagnetic ordering perpen- 
dicular to the chains, due to the fairly strong antifer- 
romagnetic couplings J ac i and J aC 2- The frustrated in- 
trachain couplings J c \ and J C 2, tending to compensate 
each other, play only a minor role, connecting next- 
nearest and more distant spins on the square lattice. 
The transition from the antiferromagnetic to the para- 
magnetic phase occurs at about Tn = 11.5±0.5K (shift- 
ing to somewhat lower temperatures when assigning the 
uniaxial anisotropy mainly to the intrachain couplings). 
This apparent agreement with the experimentally ob- 
served transition temperature of 10. 5K (Ref. [l]) resp. 
9.7K (Ref. 0) is, however, misleading. Quantum fluctua- 
tions are known to reduce the transition temperature in 
closely related two-dimensional anisotropic spin-1/2 an- 
tiferromagnetic quantum Heisenberg models by a factor 
of about 1.5 in the extreme Ising limit, as compared to 
T)v in the corresponding classical models with unit vec- 
tors; for weaker anisotropy the factor tends to increase!^ 
Therefore we conclude that the transition temperature of 
the anisotropic Heisenberg model 1 is expected to be too 
low compared to the measured Tjy. 

(ii) Well above the transition temperature, say, at 
1.5 T/v < T < 3.0 Tw, the ratio r of the static sus- 
ceptibilities (defined by the magnetizations divided by 
the magnetic fields) along the b and c axes, resp., is 
found to be almost independent of temperature, both 
in experiments^ and in the simulations, see Fig. ^ F° r 
the anisotropic Heisenberg model, we find that the ra- 
tio is also largely independent of the strength of the 
field (considering fields smaller than the spin-flop field 
at zero temperature, Hf— 1.81meV), with r « 1.02 (ap- 
proaching one at higher temperatures). The measured 
ratio is appreciably larger, r w 1.3, see Fig. 1 in Ref. 0- 
This discrepancy may be, however, misleading because 
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FIG. 1: Total magnetization M(h >c \ divided by the strength of 
the differently oriented fields Hn,,c) as a function of tempera- 
ture as obtained from Monte Carlo simulations of the classical 
version of the Hamiltonian by Matsuda et al. with a single 
ion anisotropy, Ref. EL for a system with 60 x 60 spins. Note 
that Hq,, c ) is expressed in units of meV. 
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FIG. 2: Susceptibility \b = 8M b /dH b at k B T = 0.5meV and 
0.7meV as determined in simulations of the classical version 
of the Hamiltonian by Matsuda et al. for a system with 60 x 
60 spins. 



the ratio r obtained in the simulations has to be multi- 
plied by the ratio of the squares of the corresponding 
g-factors, (ffb/gc) 2 , for a correct comparison with the 
experiments^^ Indeed, the ^-factor is anisotropic, with 
9b/gc ~ 1.14 for this material. 5 - 7 Thence the simulational 
findings for the ratio r of the magnetizations in differently 
oriented fields on the anisotropic Heisenberg model with 
a rather weak anisotropy seems to be consistent with the 
pertinent experiments. 

(iii) Applying the field along the b axis, we evaluated 
the susceptibility Xb — dMb/dHb at temperatures below 
TV. Typically, at constant temperature, xt displays as 
a function of Hf, a strong, delta-like peak at Hf, sig- 
nalling the transition from the antiferromagnetic to the 
spin-flop phase and, at a larger field, -ff^ m , an additional 
weak maximum indicating the transition from the spin- 
flop phase to the paramagnetic phase, see Fig. [21 The 
spin-flop field TT| f is nearly independent of temperature 
below about 0.7 TV, decreasing then rather sharply to 
zero as T approaches TV. The upper transition field, 
H^ m , decreases quite strongly with increasing temper- 
ature already at low temperatures. These findings on 
the anisotropic Heisenberg model are in marked contrast 
with the observed behavior^ 5 , of Xb for LasCagC^^zii. 
Experimentally, one also obtains two anomalies in Xb, but 
a sharp, cusp-like singularity with a rather low maximal 
height at a small field followed by a broad and higher 
maximum at larger fields. The location of the upper 
characteristic field depends, below TV, only weakly on 



temperature. At the lower characteristic field the break- 
down of antiferromagnetic order is observed in neutron 
scattering experiments^ this field goes to zero as one ap- 
proaches TV- Experimentally, there seems to be no evi- 
dence for a transition from the antiferromagnetic phase 
to a spin-flop phase at the lower characteristic field. The 
origin and nature of that transition has still to be clar- 
ified. At any rate, the results on Xb for the anisotropic 
Heisenberg model deviate clearly from the measured be- 
havior. 



From the comparison between the thermodynamic 
measurements^ 5 - on LasCagC^-iC^i and the simula- 
tional data on the classical version of the model sug- 
gested in Ref. Q, it follows that the model, describing 
nicely the dispersion relations, seems to be incomplete, 
missing important ingredients of the real material. Per- 
haps, as already indicated beforeji the holes induced by 
the Ca ions may play an important role. A first, sim- 
plified modelling of that aspect in the framework of an 
Ising model has been proposed recently^S Furthermore, 
the couplings of spins in neighboring ac-planes may be 
relevant. Of course, extending the model of Matsuda et 
al^ to three dimensions is not expected to resolve the dis- 
crepancy about the presence of the spin-flop phase. Fi- 
nally, structural distortions in LasCagC^-iC^i may even 
call for a description going beyond pure spin models. 
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